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ABSTRACT 

The inverse Compton (IC) scattering of ultrarelativistic electrons accelerated at the 
pulsar wind termination shock is generally believed to be responsible for TeV gamma- 
ray signal recently reported from the binary system PSR B1259-63/SS2883. While this 
process can explain the energy spectrum of the observed TeV emission, the gamma-ray 
fluxes detected by HESS at different epochs do not agree with the published theoretical 
predictions of the TeV lightcurve. The main objective of this paper is to show that 
the HESS results can be explained, under certain reasonable assumptions concern- 
ing the cooling of relativistic electrons, by inverse Compton scenarios of gamma-ray 
production in PSR B1259-63. In this paper we study evolution of the energy spec- 
tra of relativistic electrons under different assumptions about the acceleration and 
energy-loss rates of electrons, and the impact of these processes on the lightcurve of 
IC gamma-rays. We demonstrate that the observed TeV lightcurve can be explained 
(i) by adiabatic losses which dominate over the entire trajectory of the pulsar with 
a significant increase towards the periastron, or (ii) by the "early" (sub- TeV) cutoffs 
in the energy spectra of electrons due to the enhanced rate of Compton losses close 
to the periastron. The first four data points obtained just after periastron comprise 
an exception - possibly due to interaction with the Be star disk, which introduces 
additional physics not included in the presented model. The calculated spectral and 
temporal characteristics of the TeV radiation provide conclusive tests to distinguish 
between these two working hypotheses. The Compton deceleration of the electron- 
positron pulsar wind contributes to the decrease of the nonthermal power released in 
the accelerated electrons after the wind termination, and thus to the reduction of the 
IC and synchrotron components of radiation close to the periastron. Although this ef- 
fect alone cannot explain the observed TeV and X-ray lightcurves, the Comptonization 
of the cold ultrarelativistic wind leads to the formation of gamma-radiation with a spe- 
cific line-type energy spectrum. While the HESS data already constrain the Lorentz 
factor of the wind, T $S 10 6 (for the most likely orbit inclination angle i — 35°, and 
assuming an isotropic pulsar wind), future observations of this object with GLAST 
should allow a deep probe of the wind Lorentz factor in the range between 10 4 and 
I0 6 . 

Key words: acceleration of particles - binaries: gamma-rays:individual:PSR B1259- 
63 



1 INTRODUCTION 

PSR B1259-63/SS 2883 - a binary system consisting 
of a 47ms pul s ar or biting around a luminous Be star 
(| Johnston et all I1992T ) - is a unique high energy labora- 
tory for the study of nonthermal processes related to the 
ultrarelativistic pulsar winds. X-ray and gamma-ray emis- 
sion components are expected from this object due to 
the radiative (synchrotron and inverse Compton) cooling 



E-mail:dmitry.khangulyan@mpi-hd. mpg.de 



of relativis tic electrons accelerated by the wind termina- 
tion shock l|Tavani fc Aronsil 19971 ; iKirk et al.lll999l '). Gener- 
ally, the particle acceleration in this complex system can 
be treated as a scaled-down in space and time (" com- 
pact and fast" ) realization of the current paradigm of Pul- 
sar Wind Nebulae (PWN) which suggests that the inter- 
action of the ultrarelativistic pulsar wind with surround- 
ing medium leads to the format i on of a relativistic stand- 
ing shock dRees fc Gunnl 1 1974 iKennel fc Coronitil Il984l ; 
lHarding fc Gaisserlll990l ; lTavani fc Aronslll997h . In the case 
of strong and young pulsars, the shock-accelerated multi- 
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TeV electrons should give rise to observable X-ray (syn- 
chrotron) and TeV (inverse Compton) nebulae with typi- 
cal linear size ~ 0.1 — 10 pc. The unambiguous association 
of some of the recently discovered extended TeV gamma- 
ray sources with several distinct synchrotron X-ray PWNs 
generally supports this scenario of formation of nonthermal 
nebulae around the pulsars. 

In the binary system PSR B1259-63/SS 2883 one should 
expect a similar mechanism of conversion of the major frac- 
tion of the rotational energy of the pulsar to ultrarelativis- 
tic electrons through formation and termination of the cold 
electron-positron wind. On the other hand, in such systems 
the magnetohydrodynamic (MHD), acceleration and radia- 
tion processes proceed under essentially different conditions 
compared to the PWN around isolated pulsars. In particu- 
lar, due to the high pressure of the ambient medium caused 
by the outflow from the companion star, the pulsar wind 
terminates quite close to the pulsar, R ^ 10 12 cm. Con- 
sequently, in such systems particle acceleration occurs at 
presence of much stronger magnetic field (B ~ 0.1 — 1 G), 
and under illumination of intense radiation from the normal 
star with a density 
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where L s tar is the luminosity of the Be star and R is the 
distance between the acceleration site and the Be star. The 
discussed ranges of the temperature and the luminosity of 
the star SS2883 vary within T ~ 2.3 ■ 10 4 - 2.7 ■ 10 4 K and 
L star =± 3.3 ■ 10 37 - 2.2 • 10 38 erg/s. This implies that both 
the acceleration and radiative cooling timescales of TeV elec- 
trons are of order of hours, i.e. comparable or shorter than 
the typical dynamical timescales characterizing the system. 
This allows a unique "on-line watch" of the extremely com- 
plex MHD processes of creation and termination of the ul- 
trarelativistic pulsar wind and the subsequent particle ac- 
celeration, through the study of spectral and temporal char- 
acteristics of high energy gamma-radiation of the system. 
The discovery of TeV gamma-radia tion from PSR B1259 - 
63/SS2883 by HESS collaboration l| Aharonian et all 120051 ) 
provides the first unambiguous evidence of particle acceler- 
ation in such systems to TeV energies. 

Remarkably, in spite of the complexity of the binary 
system PSR B1259-63/SS2883, one may calculate with quite 
high precision the spectral and temporal features of gamma- 
radiation based on only a few model assumptions concern- 
ing, in particular, the magnetic field, the acceleration rate 
as well as the non-radiative (adiabatic or escape) losses of 
electrons. The observational uncertainties involved in these 
calculations are mainly related to the luminosity of the op- 
tical companion SS 2883. The basic free parameters used 
in calculations are the rates of particle acceleration and the 
nonradiative losses caused by adiabatic expansion and es- 
cape of electrons. 

The orbital elements of the system are well known. 
The orbit is quite eccentric (e = 0.87) with the minimum 
and maximum distances between the pulsar and the star 
Do = 9.6 x 10 12 cm and D a = 1.4 x 10 14 cm at the periastron 
and the apastron, respective ly; the inclination of orbit is 
i ~ 35° ijjohnston et alJll992h . The separation between the 
pulsar and the optical companion versus the epoch is shown 
in Fig[T] In the same figure we show the possible locations 



of pulsar passage thro ug h the stellar disk as discu ssed by 
I Johnston et"afl <|2005l ') and lChernvakova et all (|2006l ). In the 
first paper the disk location is determined by the time of dis- 
appearance of the pulsed radio emission which is explained 
by absorption of the pulsar emission in the stellar disk. How- 
ever, as long as the physics of interaction of the pulsar winds 
with the stellar disk is not firmly esta blished, alternative 
model s are not excluded. For example, IChernvakova et "all 
(2006) noticed that the maxima of lightcurves of nonpulsed 
radio, X-ray and TeV gamma-ray are quite close to each 
other, and proposed that the increase of the nonthermal en- 
ergy release happens when the pulsar crosses the disk. If so, 
the ga mma-ray emission could be result of hadronic i nterac- 
tions (|Kawachi et al J 12004 IChernvakova et al.ll2006l ). This 
hypothesis implies, however, a different location of the disk 
compared to t he one derived from the eclipse of the p ulsed 
radio emission ({Johnston et al . 2005; B ogomazovll2005l ) . and 
therefore requires an independent confirmation based on a 
stronger evidence of correlation of radio, X-ray and TeV 
fluxes, as well as detailed study of the reasons of such cor- 
relation. 

Meanwhile, the inverse Compton scattering remains 
the most plaus i ble gamma-ray produ c tion mechanism 
, m JTavani fc Aronsl Il9971 ; iKirk et all Il999l ; lAharonian et all 
l2005h . In this paper we present detailed numerical stud- 
ies of the the spectral and temporal characteristics of TeV 
gamma-ray emission within the framework of the IC model 
of TeV gamma-rays. In this context we adopt the position 
of the dis k as it is derived from the eclipse of pulsed radio 
emission ijjohnston et al.ll2005h . 

The position of the shock wave is determined by inter- 
action of the pulsar wind with the stellar wind, therefore the 
distance to the shock is a function of time. For the magnetic 
field lines frozen into the pulsar w ind, one has B o c r~^. 
It is also expected that r s h oc D (|Kirk et al.lll999l ). thus 
icb oc D~ 2 . As long as the mass flux density from a Be star 
is e xpected to be signifi cantly higher than one from the pul- 
sar (| Waters et al.|[l988h . we adopt r s h -C D, and the target 
photon density at the site of electron acceleration and radi- 
ation Wph = Lstar /(47rc)D~ 2 . Thus the synchrotron and IC 
radiation timescales have similar dependencies on the sepa- 
ration distance D, namely t Byn ,tic oc D 2 . 

The magnetic field strength in the pulsar wind at dis- 
tance r from the pulsar can be estimated as the following 
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where L s d is spindown luminosity of the pulsar; and o is 
the ratio between the Poynting and kinetic energy flux in 
the pulsar wind. While the spindown luminosity is entirely 
known and is measured to be L s d = 8.3 • 10 35 erg/s; the er 
parameter is not constrained either theoretically or observa- 
tionally. In what follows we assume the magnetic field at the 
termination shock to be B ~ 0.1 — 1G around the periastron 
epoch. This magnetic field strength corresponds to a value of 
a = 0.04 — 4 • 10~ 4 (for the distance between the pulsar and 
the termination shock r s h = 10 12 cm), which is consisten t 
with the value a = 0.02 adopted bv lTavani fc Aronsl l|l997h . 
As long as the a-parameter is rather small, we do not study 
some minor effects, e.g. magnetic field amplification on the 
termination shock, or possible impact of magnetic field ad- 
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justment on the particle flux and the Lorentz factor of the 
pulsar wind. 

For the expected magnetic field strength, the corre- 



sponding energy density is B /8tt 



Kl- 



lO 1 erg/cm 3 



The energy density of the photon field significantly exceeds 
this value (see Eq. {T}). Thus the radiation is formed in 
an environment dominated by radiation. Since the temper- 
ature of the starlight is about 2 eV, the inverse Comp- 
ton scattering proceeds in a regime with distinct features 
related to the transition from the Thomson to Klein- 
Nishina limits, depending on the scattering angle (i.e. lo- 
cation of the pulsar in the orbit ) and the electron energy 
ijKhangulvan fc Aharoniarj|2005l ). 

The TeV gam ma-ray lightcurve of b inary PSR B1259- 
63/SS2883 shows {Ah aroman et al.ll2005l ) a tendency for a 
minimum flux at the epoch close to the periastron passage, 
as well as a maximum observed 20 days after the periastron. 
Although the available TeV data do not allow robust conclu- 
sions about the lightcurve before the periastron, there is evi- 
dence of a time variable flux which indicates the existence of 
the second maximum >18 days before the periastron. While 
the light curve reported by HESS needs independent con- 
firmation by future measurements, throughout this paper 
we assume that the TeV flux decreases to wards the perias- 
tron, as stated by the HESS collaboration l|Aharonian et abl 
120051 ). Int erestingly, the X-ray observations show a s imilar 
behavior (|Tavani et al.|[l996l ; IChernvakova et al.|[2006l ). 

Below we discuss 3 different possible scenarios which 
could explain the drop of the TeV gamma-ray luminosity 
close to periastron: (i) nonradiative (adiabatic or escape) 
losses of electrons; (ii) "early" (sub- TeV) cutoffs in the en- 
ergy spectra of shock-accelerated electrons due to the in- 
crease of the rate of Compton losses, and (iii) decrease of 
the kinetic energy of the pulsar wind before its termination 
due to the Comptonization of electrons in the cold ultra- 
relativistic wind. Note that generally in the X-ray binaries 
with luminous companion stars the photon-photon absorp- 
tion may have a strong impact in the formation of the TeV 
gamma-ray lightcurves. However, in the case of PSR B1259- 
63/SS2883 the absorption effec t appears to be not significant 
l|Kirk et all 1 19991 ; |Pubusll2005l ). 



2 THE ELECTRON DISTRIBUTION 
FUNCTION 

Formally, the radiation seen by an observer is contributed 
by electrons of different ages, i.e. by electrons from differ- 
ent locations of the pulsar during its orbiting. However, 
since the radiative cooling time of high energy electrons 
(E ^ 100 GeV) is quite short (see below), we effectively 
see the radiation components from a localized part of the 
orbit with homogeneous physical conditions. This allows us 
to reduce the treatment of time-evolution of energy dis- 
tribution of electrons to the w ell-known equation (see e.g. 
iGinzburg fc Svrovatskiil l|l964l )) 
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where 7 = 7 ic + 7 sync h + 7ad;7ic, Tsynch, 7ad are electron 
energy loss rates (IC, synchrotron and adiabatic, respec- 
tively) and Q(t, 7) is the acceleration rate. The applicability 



of the Eq.((3]) to the case of IC losses in the Klein-Nishin a 
regime was shown e.g. in iKhangulvan fc Aharonianl (|2005h . 
The solution of this equation has the following form 

Toff 

n(t, 7) = ^ / Q(t - r, 7 ') e - T <™'>/ T ™dy , (4) 
7 

where Q(t, 7) is the acceleration rate at the given epoch; and 
7 c ff is implicitly defined by the following equation: 
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And for r(-y, 7) one has 
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Since the cooling time of electrons is much shorter than 
the characteristic dynamic times of the system, we can use 
the steady-state distribution function of electrons at given 
epoch t (or at given position of the pulsar in the orbit), 



Q(t,7')e- r(7 ' 7 ' )/Tosc d7', 



(7) 



where 7 max is the maximum Lorentz factor of injected par- 
ticles. 

It should be noted that at low energies this solution 
may have a limited applicability because of long radiative 
cooling time of electrons. The minimum energy of electrons 
for which the solution Eq.© remains correct, is determined 
from the condition for the cooling time: t C ooiin g < tdyn where 
tdyn ~ 1 day (the time during which the distance between 
the pulsar and the optical star, and other principal param- 
eters change less than 10%, even at the epochs close to the 
periastron). This condition gives 
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While very high energy electrons "die" due to radiative 
losses inside the acceleration region and cannot effectively 
escape, low energy electrons can escape from the source. The 
maximum energy of electrons which escape from the source 
is determined by the characteristic escape time, 
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GeV. 
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These electrons form a quasi-stationary halo around the 
binary system which can contribute to the overall IC radia- 
tion of the source. This component is to be formed in the 
Thomson regime {E 1 ~ e p h7c = (3&r)7e) arm taking into 
account Eq.([9]) one obtains a typical energy of the halo ra- 
diation, 



E 
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(Jesc 
To 3 ! 



GeV. 



(10) 



This radiation can be detected by GLAST as a quies- 
cent component. Note however that it can be significantly 
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suppressed in the case of a low-energy cutoff in the accel- 
eration spectrum of electrons as is often assumed for PWN 
in general, and fo r this binary system, in particular (see e.g. 
iKirk et all (|l999l )). 

Below we assume a power-law distribution for the ac- 
celerated electrons with an exponential high energy cutoff, 

Q(t,-/) = A-y~ a exp [-7mc 2 /£ c ,max] , (11) 

where A is the normalization coefficient related to the 
fraction of pulsar wind particles turned into a high-energy 
power-law distribution. The cutoff energy i? e ,max is deter- 
mined from the balance between the acceleration and energy 
loss rates, therefore it is a function of time. Generally, the 
total acceleration power determined by the parameter A is 
also time-dependent. 

If the radiative cooling times of electrons are smaller 
than the adiabatic and escape losses, the radiation of elec- 
trons proceeds in the saturation regime, i.e. the source works 
as a calorimeter. In this regime the energy of accelerated 
electrons is radiated away due to synchrotron and IC pro- 
cesses: 

L Byn = — J ,_l -^injection , (12) 

Isjro + tjc 

Lie = — 1 ! -^injection / ' (0) . (13) 

tsyu + t\Q 

The function f(0) is determined by the angular anisotropy 
of IC radiation, where 6 is the angle between the line of 
sight and the line connecting the pulsar and the Be star. In 
principle, the function / can vary significantly, but in the 
case of PSR B1259-63/SS2883 the variation of this function 
does not exceed two. 

The lightcurve of 1 TeV gamma-rays produced by elec- 
trons with injection spectrum given by Eq.(QTj with time- 
independent parameters A and -E e , m ax is shown in Figf2]by 
dash-dotted line. It is important to note that the almost 
time- independence of this curve, with a weak maximum just 
before the periastron, is the result of the anisotropy of IC 
scattering. Such a lightcurve which implies almost constant 
flux (within a factor of two) over the entire orbital period, 
is in obvious conflict with the HESS observations which 
show noticeable reduction of the flux towards the perias- 
tron l|Aharonian et af]|2005l ). To achieve such a behavior of 
the lightcurve, assuming a constant electron injection, one 
should introduce additional energy losses. This cannot be 
achieved by increasing the magnetic field, because it would 
lead to an increase of the synchrotron flux close to the pe- 
riastron, in contrast to X-ray o bservations (|Tavani fe Arona 
ll997l ; lc"hernvakova et al] 2006). This implies that one needs 
to introduce additional "invisible", i.e. nonradiative energy 
losses. This case is discussed in Section [3] Alternatively, one 
may reduce the flux both of synchrotron X-rays and of IC 
gamma-rays assuming a tendency of decrease of cutoff en- 
ergy in the spectrum of accelerated electrons -E c ,max, or as- 
suming reduction of the total power of accelerated electrons 
close to the periastron. These two cases are discussed in 
Sections [4] and O respectively. 
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Figure 1. Separation between the pulsar and the star versus 
time (0 corresponds to the periastron passage). The vertical dot- 
ted lines correspond to the location of interaction of the pulsar 
wind with the stellar disk bas ed on the observation s of the eclipse 
of the pulsed radio emission l| Johnston et al.ll2005h . The vertical 
dashed-dotted lines correspond to the disk location suggested by 
Chernyakova et al. (2006). 
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Figure 2. Inverse Compton gamma-ray lightcurves calculated for 
1 TeV photons. The curves correspond to three different models. 
(i) The dashed-dotted line corresponds to the constant injection 
power and energy cutoff in the spectrum of electrons, (ii) The 
dashed line corresponds to the case similar to (i), but assuming 
time-dependent energy cutoff (see Sec[4]l). The parameter n which 
characterizes acceleration rate is fixed at the level r\ = 7 ■ 10 3 ; the 
magnetic field varies with the separation distance as B(D) ex 1/D 
with Bo = 0.05 G at periastron. (Hi) The solid line corresponds 
to the case similar to (ii) with an additional assumption about 
escape of particles with escape time T CS c = 5 ■ fO 3 s). The exper- 
imental points are from the HESS observations at energy 1 TeV 
dAharonian et alj|2005h . 

3 NONRADIATIVE LOSSES 

The interaction of the pulsar wind with the ambient medium 
results in a complex shock wave structure where adiabatic 
and escape losses may play a dominant role. Indeed, the 
escape time T osc can be as short as 

D J 10 3 s (for periastron) , . 

csc = {cJS) ~ \ 10 4 s(forapastron) ' ! " ' ' 

i.e. comparable or shorter (depending on energy) than the 
radiative cooling time of electrons. The adiabatic losses 
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(t ~ R/v, where R is the characteristic radius and v is the 
speed of expansion of the the emission region) can be even 
faster since R <C D and the source can expand relativisti- 
cally. We note, that although under certain conditions the 
particle escape and the adiabatical loss times may be rather 
connected, in general these two timescales are indeed dif- 
ferent. For example, in the case of fast diffusion, electrons 
diffuse away from the source, without suffering any adiabat- 
ical losses. 

Both adiabatic and escape timescales cannot be calcu- 
lated from the first principles given the complexity of the 
systerrQ. Instead, the characteristic timescales of nonradia- 
tive losses can be derived phenomenologically, namely from 
the observed gamma-ray lightcurve. As discussed above, 
in order to explain the reduction of both X-ray and TeV 
gamma-ray fluxes close to the periastron, one should in- 
crease the rate of nonradiative losses. The reported TeV 
gamma-ray observations are quite sparse and unfortunately 
allow a broad range of lightcurves. In Fig. [3] we show an 
example of a lightcurve which matches the HESS data. The 
time-profile of the rate of nonradiative losses derived from 
the lightcurve along with the additional assumption that 
10% of pulsar spindown luminosity is converted (through 
the termination shock) into relativistic electrons with con- 
stant (for all epochs) injection rate and energy spectrum 
given by Eq. (|ll[) with -E e ,max = 10 TeV, a — 2 is shown 
in Fig(3] (small panel). For magnetic fields we assume a 
B(r) (x D^ 1 dependence with B = 0.1 G at periastron. 
To reconstruct the adiabatic loss profile, we first performed 
calculations with the "best guess" initial profile for adia- 
batic losses. Namely, if we assume that the TeV lightcurve 
has a minimum around periastron, then obviously the adi- 
abatic losses should increase closer to the periastron. Also, 
in order to avoid an overestimate of TeV fluxes compared to 
the fluxes (or upper limits) reported by HESS at phases well 
beyond periastron, one should assume some incr ease of adi- 
abatic losses at large separations (see as well in Kirk et al.l 
1120051 )). Then we "improve" the shape of this profile using 
several iterations. 

It can be seen that in order to match the lightcurve 
shown in Fig. [3] one needs a very sharp increase of adiabatic 
losses with characteristic time T ~ 100 sec. This can be nat- 
urally related to a much smaller size of the emission region at 
periastron, i.e. the region occupied by relativistic electrons 
accelerated by the termination shock is a denser region closer 
to the star. In the case of termination of the wind in the 
highest density environment which coincides with the pas- 
sage of the pulsar through the stellar disk we would expect 
even higher nonradiative losses, so one should expect some 
deviation from the adopted smooth symmetric profile of adi- 
abatic losses shown in Fig(3] (the small panel) . Interestingly, 
the TeV gamma-ray flux at t — 20 days from periastron is 
lower in comparison with the adopted reference lightcurve. 
Note, that this "anomaly", related to the four points after 
periastron passage, appears also in other models (see be- 
low). This "anomaly" can be interpreted as a result of the 
enhanced nonradiative losses in the disk. Another natural 



Concerning this complicated issue, we would like to note a 
possibility to get this in formation by a direct numerical simulation 
iBoeovalov et al.ll2007l) . 
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Figure 3. Main panel: The lightcurve of 1 TeV gamma-rays de- 
tecte d by HESS from PSR B1259-63/SS2883 jAharonian et al.l 
120051) . A reference lightcurve adopted for derivation of the time 
profile of nonradiative energy losses of electrons is also shown. The 
two vertical gray zones correspond to the position of the stellar 
disk. The somewhat lower flux of gamma-rays at t ~ 15 days after 
the periastron can be associated with the enhanced losses in the 
disk, and thus may cause more irregular profile of energy losses. 
Small panel: The reconstructed time profile of adiabatic energy 
loss rate derived for the reference TeV lightcurve shown in the 
main panel (see the text). 

reason for the reduction of gamma-ray emissivity could be 
the deficit of the target photons in the stellar dislfl How- 
ever, the large statistical and systematic errors of TeV fluxes 
do not allow certain conclusions in this regard. Therefore the 
reference lightcurve in Fig. [3] should be treated as a reason- 
able approximation for derivation of basic parameters of the 
system. 

At epochs far from periastron the rate of the required 
nonradiative losses drops significantly, however it still re- 
mains faster than IC and synchrotron losses with character- 
istic time T ~ 10 4 sec (see FigfJJ). Note that the curves in 
FigS] correspond to electrons of energy 1 TeV. However, as 
it is seen from FigO nonradiative losses should dominate 
over radiative losses at all energies of electrons and during 
the entire orbit of the pulsar. 

In Fig[6]we show lightcurves calculated for four differ- 
ent gamma-ray energies. The lightcurve for 1 TeV coincides, 
by definition, with the reference lightcurve shown in Fig(3] 
The lightcurves generally are similar what is explained by 
a dominance of adiabatic losses at all electron energies. At 
the same time the energy spectra of IC radiation at dif- 
ferent energy bands are significantly different. Indeed, the 
dominance of adiabatic or energy-independent escape losses 
maintains the acceleration spectrum of electrons unchanged. 
Thus at energies E -C Smax, the gamma-rays produced in 
the Thompson regime (E n ^ 10 GeV) will have a power-law 
spectrum with photon index (a + l)/2, while in the deep 

2 Even assuming that the photons from the star absorbed in the 
disk are fully re-radiated at other wavelengths, one should never- 
theless expect significant reduction of the photon energy density 
in the disk because of the isotropisation of the initial radial dis- 
tribution of photons from the star. 
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Figure 4. The energy losses for 1 TeV electrons versus the sep- 
aration distance between the pulsar and the companion star. 
The solid line is for synchrotron losses for B = 0.1G at peri- 
astron, and B(D) oc D _1 ; the dashed-dotted line is for IC losses 
for black body distribution of target photons with temperature 
T = 2.3 ■ 10 4 K diluted with coefficient n = (R ic /2D) 2 , where 
is the radius of the star. The dashed line corresponds to adiabatic 
losses for the reconstructed time-profile shown in Fig|3](the small 
panel). 
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Figure 5. The energy loss rates calculated for two epochs: at 
periastron (top panel) and at t = 35 days from periastron (bottom 
panel). The assumed parameters are the same as in Fig|4] 



Klein-Nishina regime the spectrum will be proportional to 
E~ (a+1) In E-,. This effect is seen in Fig[7] where we show the 
broadband spectral energy distribution (SED) of radiation 
at different epochs, consisting of synchrotron and IC compo- 
nents. In Fig. [8] we also show the gamma-ray spectra aver- 
aged over three periods of HESS observations in February, 
March and April 2004. Within the statistical and system- 
atic uncertainties, the a greement with the flux es reported 
by HESS is satisfactory (lAharonian et al.ll2005l ). 

For the adopted key parameters, in particular for the 
electron adiabatic loss rates shown in Figs. [4] [5] the gamma- 
ray fluxes are not sensitive to the magnetic field strength as 
long as it does not exceed 0.5 G. At the same time, the 
synchrotron flux is very sensitive to the B-field (oc B 2 ). 
It is seen in Fig[7] that for the assumed magnetic field 
0.1G the fluxes of synchrotron X-rays are significantly be- 
low 1 0~ n erg/cm 2 s, thus the y cannot explain the X-ray 
data (|Chernvakova et al.ll20od ). On the other hand, assum- 
ing a specific value of magnetic field £?» = 0.45G at peri- 
astron, one can increase the synchrotron X-ray flux to the 
observed level, and, at the same time keep the gamma-ray 
fluxes practically unchanged. Not that any significant (20 
% or so) deviation from this value of the magnetic field 
would lead to the reduction of X-ray fluxes (for B ^ B*), 
or gamma-ray fluxes (for B ^ B*) below the observed flux 
levels. The X-ray lightcurve, calculated for the B-field de- 
pendence B = 0.45(D /£>)G, i s shown Figg] While the cal- 
culated lightcurve is in reasonable agreement with observa- 
tions around the periastron and at high separations (more 
than 200 days), it significantly exceeds the fluxes between 
200 days to several days before periastron. This implies a 
much weaker magnetic field, and thus the synchrotron radi- 
ation fails to explain the X-ray emission. A more speculative 
explanation could be that, the energy release in nonthermal 
particles for this period is suppressed, e.g. due to the the in- 
teraction with the stellar outflow. If the X-ray and gamma- 
ray components are produced by the same electron popula- 
tion, this effect will have a similar impact on the gamma- 
ray lightcurve, namely it should suppress significantly the 
gamma-ray fluxes. However the lack of gamma-ray data for 
this period does not allow any conclusion in this regard. 



4 MAXIMUM ENERGY OF ELECTRONS 

It is convenient to present the acceleration time of electrons 
in the following form: 

c 



= « 0.11 £ T ov B^vs, 



(15) 



where ?~l is the Larmor radius, and rj is a dimensionless 
constant; r\ — 1 corresponds to the maximum possible rate 
of acceleration allowed by classical electrodynamics. It is 
well known theoretically, that in case of nonrelativistic par- 
allel shocks (with shock velocity v ), r\ > (c/v) 2 >• 1 (see 
e.g. IPesses et all l| 19821 ): Ijokipiil (|2004r i: lAharonian et all 
(2002)). Although this parameter is not consistently con- 
strained by theory, r\ can significantly exceed 1 even in the 
case of relativistic shocks. In the case when the energy losses 
of electron are dominated by synchrotron cooling, the max- 
imum energy of the synchrotron radiation depends only on 
the rj parameter (see below for more details). Thus the mea- 
surements of the synchrotron radiation in the cutoff region 
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Figure 6. Gamma-rays lightcurves in the nonradiative loss dom- 
inated scenario are calculated for four different energies: Ej = 
ITeV (solid line), E 1 = 0.5 TeV (dotted line), E~, = 0.1 TeV 
(dash- dotted), = lOGeV (dashed line). The HESS measure- 
ments (Aharonian et al J2005h of 1 TeV gamma-ray fluxes are also 
shown. The model parameters are same as in Fig(4] 
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Figure 7. Broadband spectral energy distribution of synchrotron 
and IC components of radiation in the nonradiative loss dom- 
inated scenario. The spectra correspond to t = —100 days 
(dashed), -10 days (dash-dotted), days (solid), +20 days (dot- 
ted line) epochs (periastron is at t = 0). The model parameters 
are same as in Figf4] 



can give us quite robust information about rj. However, this 
simple relation does not work in the case of PSR B1259-63, 
at least for the models discussed here which assume that 
electrons are cooled via inverse Compton scattering or due 
to adiabatic losses. On the other hand, the spectral shape 
of fluxes of IC 7-rays, appear quite sensitive, especially at 
TeV energies, to the value of r). Thus if the suggested mod- 
els describe the gamma-ray production scenarios correctly, 
we can derive information about r\ from the comparison of 
model calculations with the observed gamma-ray spectra. 

In Fig llOl we show characteristic acceleration times for 
3 different values of 77 = 4 x 10 3 , 10 3 , 10 2 , together with 
synchrotron and Compton cooling timescales calculated for 
the epoch of the periastron assuming the magnetic field B = 
0.05 G. In Fig llOl the energy-independent escape time, which 
was assumed to be 10 4 s, is also shown. 




E,eV 

Figure 8. The time-averaged (over the HESS observation peri- 
ods) TeV gamma-ray spectrum shown together with HESS mea- 
surements. The model parameters are same as in Figf4] 
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Figure 9. The X-ray lightcurve calculated for the adiabatic 
loss time-profile shown in Fig|3] and magnetic field strength 
B = 0.45(D/D )~ 1 G . The point sets A1- A4, X1-X10, SI 
correspond to ASCA |Hiravama et aD Il996f). XMM-Newton 
llChernvakova et alj I2009T and BeppoSAX Ichernvakova et al.l 
120061 ) observations, respectively. The open point corresponds to 
the 20-80 keV hard X-ray flux, F x ~ 3 X 10 -11 erg/cm 2 s, as 
reported by the INT EGRAL team for the period 14.1-17.5 days 
after the periastron <|Shaw et al.ll2004h . 



The maximum energy of electrons is determined from 
the balance of particle acceleration and loss rates, in Fig llOl 
this energy is defined by the intersection of curves corre- 
sponding to the acceleration and loss times. Because of es- 
sentially different energy dependences of characteristic en- 
ergy loss times t ayn , tic and t eS c, the maximum electron en- 
ergy is determined, depending on the value of rj, by IC losses 
(a) or by escape (b) or by synchrotron losses (c) (see Fig llOl) . 

If the electron energy losses are dominated by syn- 
chrotron cooling in the magnetic field Bq = B/l G with 
characteristic time 



t syn w 400Bo 2 E^ v s, 

the corresponding maximum energy of electrons is 



(16) 
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Figure 10. The acceleration and cooling times of electrons at 
periastron. The solid line corresponds to the IC cooling time ob- 
tained with accurate numerical calculation; the thick solid line 
is the IC cooling time given by the Eq. 1181 : the dashed line 
corresponds to the synchrotron cooling time (B = 0.05 G); the 
dashed-dotted line is the escape time; the dotted lines are accel- 
eration times for rj = 4 ■ 10 3 (a),1J = 10 3 (b),r? = 10 2 (c). 
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Figure 11. The acceleration and energy loss times of electrons of 
different epochs. The combined IC + synchrotron cooling times 
are calculated assuming B = 0.05 (D /D) G and T = 2.3- 10 4 K. 
The acceleration times are given by Eq. 1151 assuming r} = 4 • 10 3 . 
The horizontal dash-dot-dotted line corresponds to the T eBC = 5 • 
10 4 s. The solid lines correspond to the periastron epoch; the dash- 
dotted lines to ±10 days; dashed lines to ±20days and dotted lines 
to ±100days. The maximum injection energy due to radiative 
cooling is determined by crossing of the same type lines. The 
maximum energy related to the electron escape is determined by 
crossing of an acceleration line with the escape line. 



E C ,^^2B~ 1/2 ( I ^)~ 1/2 TeV - ( 17 ) 

Note that in the case of rj — const the maximum energy of 
synchrotron photons does not depend on the magnetic field 
(-Esyn.max oc -Ee,max -B=const), but depends on rj, namely 
£7 S yn,max ~ 100(?7/10 3 ) _1 keV. This relation contains unique 
information about the acceleration rate through the rj pa- 
rameter. 
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Figure 12. Maximum energy of electrons at the pulsar wind ter- 
mination shock. The solid line corresponds to the cutoff energy 
caused by radiative cooling. The dashed line corresponds to the 
case of combined radiative and escape losses. The model param- 
eters are same as in Fig llll 

In the regime when IC cooling dominates over the syn- 
chrotron cooling, the radiation cooling time is determined 
by 

tic « 7- 10 s Wq 1 £tcv s, (18) 

where wo is the energy density of the target photons in 
erg/cm 3 units. In Fig llOl we show the accurate numerical 
calculation of the IC cooling time. It can be seen that above 
1 TeV Eq. (|18[) provides quite accurate approximation of the 
IC cooling time. 

The corresponding maximum energy of accelerated elec- 
trons is 

£ e , m ax ~ 9 • 10 5 (B G /w f s (^L) 33 TeV . (19) 

This somewhat unusual dependence of -E G , ma x on the pho- 
ton density wo is the result of IC scattering in deep Klein- 
Nishina regime. Obviously, in the Thomson regime E e max C< 
(Bq/wo) 1 r\~ x l' 2 . The very strong dependence of £ e ^ max in 
Eq. (|19|l on three highly variable parameters, B, w and rj, 
allows variation of i5 ,max in very broad limits. For exam- 
ple, for the B oc 1/D type dependence of the B-field, and 
assuming constant rj, the increase of the separation between 
the compact object and the star by a factor of two would 
lead to the change of i? e ,max by a factor of 2 3 ' 3 ~ 10, and 
correspondingly to dramatic variation of the flux of highest 
energy gamma-rays. 

Finally, the escape of electrons may also have a strong 
impact on the variation of i? e ,max depending on the position 
of the pulsar. Actually the effective escape of electrons from 
the acceleration site is somewhat shorter than the escape 
time. So in this case, ignoring the radiative energy losses 
of electrons and taking the upper limit for the escape from 
accelerator, one has 




where T is the escape time of electrons. 

Obviously, all relevant timescales depend on the pulsar 
position in the orbit, therefore the high energy cutoff in the 
spectrum of electrons is expected to be variable. As long 
as there no theoretical predictions for possible rj parame- 
ter dependence on physical conditions in the accelerator, in 
what follows we assume rj to be constant. In Fig llll we show 
the radiation and acceleration timescales for different epochs 
- at periastron and ±10,20, 100 days from the periastron. 
For the chosen model parameters, B — 0.05(Z?o/ D) G and 
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rj = 4x 10 3 , the cutoff in the eiectron spectrum at the perias- 
tron is determined by IC fosses, whiie at farge separation dis- 
tances the synchrotron and escape fosses pfay the more im- 
portant roie in formation of the cutoff. This is demonstrated 
in Fig |12l where the high energy cutoff in the eiectron spec- 
trum is shown as a function of epoch. Sofid fine corresponds 
to the case of radiation (fC and synchrotron) losses, fn this 
case one expects significant reduction of the cutoff energy at 
epochs cfose to the periastron, where strong fC fosses push 
the cutoff energy down to ^ 1 TeV. Far from the perias- 
tron, the cutoff energy can increase up to 10 TeV, uniess the 
losses due to escape become dominant. The IC coofing time 
at the epoch with separation D is t coo i ~ 10 3 s(D/Do) 2 s. 
Therefore, if the characteristic escape time is about 10 4 s, 
the impact of particie escape becomes important for separa- 
tions D ^ 3 Do- This effect is demonstrated in Fig |12l where 
(time and energy-independent) escape time T CBC = 5 x 10 4 s 
is assumed. One can see that for chosen modef parame- 
ters the cutoff energy is a weak function of time with a 
local minimum (~ 0.5 TeV) at periastron, and two maxima 
(~ 2.5 TeV) at ±20 days. 

It is important to note that the introduction of es- 
cape fosses is crucial for expfanation of the observed TeV 
lightcurve in this scenario. Indeed, while the reduction of 
the cutoff energy in the spectrum of eiectrons due to en- 
hanced IC losses satisfactoriiy expfains the minimum at the 
periastron, this wouid impiy much higher fluxes at large 
separations in contrast to the HESS observations. The ad- 
ditional assumption that electrons suffer afso significant es- 
cape losses (Tcsc = 5x 10 4 s) allows dramatic suppression of 
the gamma-ray fluxes beyond \t\ > 20 days (compare dashed 
and dot-dashed curves in Fig(2|. 

The impact of the variation of refative contributions of 
radiative and escape fosses on the formation of the energy 
distributions of efectrons is demonstrated in Fig |13l The cor- 
responding iightcurves of inverse Compton gamma-rays at 
E 1 = 1 TeV, 500 GeV, 100 GeV and 10 GeV, and 1-10 
keV synchrotron photons are s hown in Fig ll4l For compar- 
ison the HESS measurements (|Aharonian et al ]|2005|) of 1 
TeV gamma-ray fluxes are afso shown. The agreement of 
calculations with the HESS lightcurve is rather satisfactory 
except for somewhat higher predicted flux at the epoch of 2 
weeks after the periastron which coincides with the puisar 
passage through the steffar disk. In Fig |15l we compare the 
energy spectrum reported by HESS with the average TeV 
gamma-ray spectrum calcufated for the period of the HESS 
observations in February 2004. Aithough it is possibie to 
achieve a better agreement with the measurements, at this 
stage the attempt for a better spectrai fit coufd be hardfy 
justified given the statistical and systematic uncertainties of 
the measurements. 

Through a variation of -E e ,max the fightcurves at TeV 
and GeV energies can have quite different profiles. Namely, 
the TeV Iightcurves have a clear minimum at periastron 
which is explained by the sub- TeV cutoff in the spectrum 
of accelerated electrons. At the same time this cutoff in 
the eiectron spectrum is still sufficientiy high and therefore 
does not have a strong impact at GeV energies. Therefore 
the GeV Iightcurves show their maximum a few days before 



the periastron 0. It is important to note that the significant 
drop of gamma-ray fluxes at large separations is due to the 
escape losses, otherwise one should expect rather constant 
flux with a weak maximum close to periastron. Also this 
model predicts different energy spectra of gamma-ray below 
100 GeV at different epochs. Indeed, at large separations, 
when the escape fosses dominate, the injection spectrum of 
efectrons remains unchanged, therefore we expect noticeably 
harder gamma-ray spectra in the GeV energy band at epochs 
\t\ > 20 days (see Fig fTo) . 

Remarkabfy, the calculated fluxes at GeV energies are 
well above the sensitivity of GLAST which makes this source 
a perfect target for future observations with GLAST. It 
should be noted, however, that the fluxes at GeV energies 
couid be significantly suppressed due to a possible low en- 
ergy cutoff in the acceleration spectrum of efectrons. This 
is a standard assumption in the models of PWN, see e.g. 
iKennet fc Coronitil i| 19841 ); iKirk et al l (|l999h - 

In this scenario, the magnetic field energy density at 
the shock shouid be significantfy less than the energy den- 
sity of steflar photons. Thus, assuming the same strength of 
the magnetic fieid in the acceleration and radiation regions, 
one obtains quite low synchrotron fluxes (see Fig. I14|l . This 
would imply that the obser ved X-rays have non-synch rotron 
origin (e.g. IC origin; see (|Chernvakova et al.l [20 06) ) . An- 
other possibility is to assume that the magnetic field in the 
radiation region is somewhat higher than at the shock (note 
that that a similar situat ion takes place in the Crab nebula 
|Kennel fc Coroniti| [l984'l'l. In Fig[T7]we show the lightcurve 
of 1-10 keV X-rays, assuming that in the radiation region 
the magnetic field is stronger by a factor of 8. If so, the syn- 
chrotron X-ray flux could achieve the observed flux level. In 
this scenario the X-ray and gamma-ray production regions 
are essentially different, although they could partly overlap. 
While the bulk of X-rays is formed in a magnetized region(s) 
far from the shock (where the electrons are accelerated) , the 
IC gamma-rays come from more extended zones, which in- 
clude also the site of particie acceleration. In Fig |17l we show 
also the 20-80 keV hard X-ray flux as reported by the IN- 
TEGRAL team fo r the 14.1-17.5 days after the periastron 
IShaw et a.1.1 (|2004l ) . Note that these days correspond to the 
minimum of the gamma-ray fluxes which we explain quan- 
titatively by the passage of the pulsar through the disk, e.g. 
due to the the deficit of the target radiation field. 



5 COMPTONIZATION OF THE UNSHOCKED 
WIND 

While in the previous sections we tried to explain the ob- 
served modulation of TeV gamma-rays fluxes by energy 
losses of acceferated electrons at the termination of the wind, 
it is interesting to investigate whether one can relate the ob- 
served TeV lightcurve to the Compton losses of the kinetic 
energy of the bulk motion of the cold ultrarelativistic wind. 

Generally this effect can be important in a binary sys- 
tem with a high luminosity companion star. Although the 

3 Note that the shift of the position of the maximum is caused, 
as discussed above, by the anisotropy of the Compton scattering, 
but not by the change of the target photon density as long as the 
IC proceeds in the "saturation regime". 
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Figure 13. Electron energy distributions at different epochs 
(0, ±10, ±20, ±100 days to periastron passage). The model pa- 
rameters are same as in Fig llll 
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Figure 15. The averaged spectrum of IC gamma-rays compeared 
with the spectrum measured by HESS during the period Feb. 
2004. The model parameters are same as in Fig |14l 
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Figure 14. The gamma-ray lightcurves expected in the scenario 
of variation of the energy cutoff. The solid line: E 1 = 1 TeV, the 
dotted line: _E 7 = 0.5 TeV, the dash-dotted line: _E 7 = 0.1 TeV, 
the dashed line: Ej = 10 GeV and the dash-dot-dot-dotted 
line: Ej = 1 — lOkeV (synchrotron). The model parameters 
are same as in Fig llll The injection spectrum was assumed to 
be Q(E e ) oc i?^~ 2 exp (— i? e /E e ,max). The injection acceleration 
rate of electrons was assumed to be at level of 5% of the pulsar 
spind own luminosity. The HESS measurements ( Aharonian et aD 
2005) of 1 TeV gamma-ray fluxes are also shown. The vertical 
shadowed zones correspond to the stellar disk location. 



electrons in cold pulsar winds do not suffer synchrotron 
losses, a significant fraction of original kinetic energy of the 
bulk motion of the electron-positron wind could be radiated 
away due to Comptonization. Thus the power available for 
acceleration of electrons at termination of the wind depends 
on the position of the pulsar. Obviously, in this scenario we 
expect a minimum flux of gamma-rays closer to periastron. 
In other words, while in the previous section the modula- 
tion of the gamma-ray flux is linked to the i? e ,max, in this 
scenario the gamma-ray flux variation depends on the pa- 
rameter A characterizing the acceleration power of electrons 
given by Eq.([TT}. 

According to th e standard PWN model 
l|Kennel fc Coronitil Il984l ) the isotropic cold electron- 
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Figure 16. Gamma-ray spectra of synchrotron and IC radiation 
at -100 (the dashed line), -10 (the dotted-dashed line),0 (the solid 
line), +20 (the dotted line) days relative to the periastron passage. 
The model parameters are the same as in Fig |14l 



positron winc0 has a typical bulk motion Lorentz factor 
F oc 10 4 — 10 6 , thus the interaction of the wind electrons 
with starlight in the Klein-Nishina limit should lead to the 
formation of a narrow gamma-ray component with typical 
energy Ym e (? . This effect also leads to the modulation of 
the bulk motion Lorentz factor as shown in Fig ll8l The cal- 
culations in Fig |18l were performed for L s tar = 2.2- 10 38 erg/s 
and different values of the initial Lorentz factor of the pulsar 
wind. Note that in the Thomson regime i coo i oc r~ , i.e. 
the decrease of wind Lorentz factor leads to the increase of 
the cooling time. On the other hand, in the Klein-Nishina 
regime the cooling time increases with Lorentz factor as 
icooi oc F ' 7 . Therefore the maximum effect is achieved in 



4 He re we follow a model suggested by iKennel fc Coronitil 
i 19841 ). assuming an isotropic pulsar wind, but it is wor- 
thy to note that the pulsar wi nd can be strongly anisotropic 
(Bogovalov & Khangoulian 2002), thus the non- typical lightcurve 
can be a result of the interaction of two anisotropic winds. 
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Figure 17. The X-ray lightcurve calculated for the scenario with 
early sub TeV cutoff. At shock front the magnetic field was as- 
sumed to be Bg = 0.03 (Dq/D) G, and in the emitting region to 
be higher by a factor of 10. The acceleration efficiency has value 
r\ = 4 X 10 3 . An orbital-dependent escape time was assumed to 
be i osc = 1.75 x 10 4 (D/D g) 1/2 s. The point sets A1-A4, X1-X10, 
SI correspond to ASCA Jffiravama et al1ll996ri. XMM-Newton 
llChernvakova et al.l l2006h and BeppoSAX l lChernvakova et all 
2006) observations, respectively. The open point corresponds to 
the 20-80 keV hard X-ray flux, F x ~ 3 X 10 -11 erg/cm 2 s, as 
reported by the INT EGRAL team for the period 14.1-17.5 days 
after the periastron (Shaw et al. 20041). 



the Thomson-to-Klein-Nishina transition region, i.e. around 
r = 10 5 . 

As it follows from Fig llBI the initial Lorentz factor of 
the wind for To = 10 is reduced by 40% at periastron. In- 
terestingly, minimum reduction of the initial Lorentz factor 
(~ 5%) happens around t = ±20 days (see in Fig l 18f) . Since 
the kinetic energy of the wind radiated away due to the 
Comptonization is determined by the starlight density and 
the length of the unshocked wind Al — r s h, the lightcurve is 
explained by the combination of two factors - dependence 
of the starlight density on the separation D, and the dis- 
tance to the wind termination point (it is assumed that the 
gas density of the stellar wind decreases as D 2 ). Obviously, 
in the case of electron-positron pulsar wind, the kinetic en- 
ergy of bulk motion of the wind, and consequently the rate 
of shock accelerated electrons have similar time behaviors: 
A — AoT(t)/To, where Ao characterizes the original power 
of the wind. 

Although qualitatively this behavior agrees with the 
TeV lightcurve detected by HESS, the effect of reduction 
of the kinetic energy of the wind is not sufficient to ex- 
plain quantitatively the observed TeV lightcurve. Indeed, 
the Comptonization of the wind can lead to the reduction 
of the energy flux of TeV gamma-rays at periastron by only 
a factor of ^ 1.5, while the HESS observations show more 
significant variation of the gamma-ray flux. Assuming some- 
what larger, by a factor of two, luminosity of the optical 
star, one can get a better agreement with the observed TeV 
lightcurve. However, the range of luminosity of the star dis- 
cussed in the literature favors a lower lu minosity of the star 
l|Tavani fc Aronslll997l ; iKirk et alJll999h . Therefore the ef- 
fect of Comptonization of the pulsar wind cannot play, even 
for To ~ 10 , a major role in the formation of the TeV 
lightcurve. 

Even so, this effect cannot be ignored in the calcula- 
tion of the overall gamma-radiation of the system. Namely 
the Comptonization of the ultrarelativistic pulsar wind un- 
avoidably leads to an additional component of gamma-rays 



produced at the pretermination stage of the wind. Due to 
the inverse Compton scattering of monoenergetic electrons 
on target photons with narrow, e.g. Planckian distribution, 
we should expect a specific, especially for Fp ^ 10 , line- 
type gamma-ray e mission l|Bogovalov fc Aharonian! |2000| ; 
IBall fc Kirkll2000l ; IBall fc DoddH200ll ). IBall fc Kirkl (|2000f ) 
have calculated the IC radiation of the freely expanding 
wind, i.e. under assumption that the wind was not termi- 
nated. These calculations can be hardly applied for this sys- 
tem, since we see broad-band IC radi ation of the shocke d 
wind. As it was shown in the paper bv lBall fc Doddl i|200ll ). 
such an assumption leads to significant (an order of magni- 
tude) overestim ation of the gamma-r ay flux. Unfortunately 
in the paper by IBall fc Doddl l|200ll ) there are only calcu- 
lations for one value of Lorentz factor T = 10 6 . To supply 
this gap we performed calculations for a number of proba- 
ble Lorentz factors. The results of calculations of gamma-ray 
spectra of the unshocked wind are shown in Fig |19l Compar- 
ison of these calculations with the average energy spectrum 
of PSR B 1259-63 measured by HESS excludes the initial 
Lorentz factor of the wind To = 10 6 . Otherwise, the flux of 
the Comptonized emission of the unshocked wind would ex- 
ceed the observable flux even at apastron (see Fig |19f) . Due 
to the energy range of gamma-rays from this source available 
for HESS (E > 300 GeV), future observations unfortunately 
cannot significantly improve this limit. On the other hand, 
such studies can be performed by GLAST the sensitivity of 
which seems to be adequate, as is shown Fig |19l for a deep 
probe of the initial wind Lorentz factors of PSR B 1259-63 
within 10 4 to 10 6 . Thus, GLAST has a unique potential to 
prove the current pulsar wind paradigm which assumes that 
the bulk of the spin-down luminosity of the pulsar is trans- 
formed to a cold wind with Lorentz factor exceeding 10 4 . 



6 SUMMARY 

One of the recent exciting results of observational 
gamma-ray astronomy is the detection of TeV gamma- 
ray signal from the bi nary system PSR B1259-63/SS2883 
{Ah aroman et al.ll2005l ). While the absolute fluxes and en- 
ergy spectra of TeV emission detected by HESS can be ex- 
plained quite_weU hithe framework of inverse Compton 
model (|Kirk et al.lll999T ). the observed TeV lightcurve ap- 
pears to be significantly different from early predictions. 
This can be considered as an argument in favor of alter- 
native (hadronic) models which relate the maximum in the 
TeV lightcurve observed after approximately three weeks of 
the periastr on to the interaction of the pulsar wind with th e 
stellar disk (jKawachi et al.ll2004l ; Ichernvakova et al.ll2006f l. 
However, the hadronic models require a rather specific loca- 
tion of the stellar disk, which needs to be confirmed, since 
this assumption does not agree with the location of the stel- 
lar disk derived fro m the pulsed radio emission observations 
l|Bogomazovll2005l ). 

The main objective of this paper is to show that the 
HESS results can be explained, under certain reasonable as- 
sumptions concerning the cooling of relativistic electrons, by 
inverse Compton scenarios of gamma-ray production in PSR 
B1259-63/SS2883. Namely, we study three different scenar- 
ios of formation of gamma-ray lightcurve in the binary sys- 
tem PSR B1259-63/SS2883 with an aim to explore whether 
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one can explain the observed TeV lightcurve by electrons ac- 
celerated at the pulsar wind termination shock. The natural 
target for the inverse Compton scattering in such a system is 
the thermal radiation from the optical star. Since the basic 
parameters characterizing the system are well known, the 
predictions of gamma-ray fluxes at different epochs can be 
reduced to a calculation of the energy distribution of the rel- 
ativistic electrons under certain assumptions concerning the 
acceleration spectrum of electrons and of both their nonra- 
diative (adiabatic and escape) and radiative (Compton and 
synchrotron) energy loss origin. In particular, we demon- 
strate that the observed TeV lightcurve can be explained (i) 
by adiabatic losses which dominate over the entire trajec- 
tory of the pulsar with a significant increase towards peri- 
astron, or (ii) by the variation of the cutoff energy in the 
acceleration spectrum of electrons due to the modulation of 
rate of inverse Compton losses depending on the position 
of the pulsar relative to the companion star. Our mod- 
els failed to explain the first four data points obtained just 
after periastron, what can be explained by a possible inter- 
action with the Be star disk, which introduces additional 
physics not included in the presented model. Although we 
deal with a very complex system, we demonstrate that the 
observed TeV lightcurve can be naturally explained by the 
inverse Compton model under certain physically well jus- 
tified assumptions. Unfortunately, the large systematic and 
statistical uncertainties, as well as the relatively narrow en- 
ergy band of the available TeV data do not allow robust 
constraints on several key model parameters like the mag- 
netic field, escape time, acceleration efficiency, etc. This also 
does not allow us to distinguish between different scenarios 
discussed above. In this regard, the future detailed obser- 
vations both in MeV/GeV and TeV bands by GLAST and 
HESS closer to periastron, as well as at the epochs when 
the pulsar crosses the stellar disk, will provide strong insight 
into the nature of this enigmatic object. Equally important 
are the detailed observations of X-rays, e.g. with Chandra, 
XMM and Suzaku telescopes. The analysis of gamma and 
X-ray data obtained simultaneously should allow extraction 
of several key parameters characterizing the binary system. 

Finally, although the Compton cooling of the unshocked 
electron-positron wind does have significant impact on the 
formation of the TeV lightcurve, the specific, line type 
gamma-radiation caused by the Comptonization of the cold 
ultrarelativistic wind should unavoidably appear either at 
GeV or TeV energies depending on the initial Lorentz fac- 
tor of the wind. Detection of this component of gamma- 
radiation, in particular by GLAST, of the unshocked wind 
will provide unique information on the formation and dy- 
namics of pulsar winds. 
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Figure 18. The time-evolution of the wind Lorentz factor just 
before the termination. The results are presented for the initial 
Lorentz factors To = 10 5 (solid line), To = 10 6 (dotted line) 
and for To = 10 4 (dashed-dotted line) for the luminosity of the 
companion star L a tar = 2.2 ■ 10 38 erg/s. Dashed line shows the 
time evolution of the wind Lorentz factor for To = 10 5 calculated 
the luminosity of the companion star -L s t a r — 3.3 ■ 10 37 erg/s. 
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Figure 19. The energy spectra of gamma-rays due to the 
Comptonization of the unshocked wind at the periastron (solid 
lines) and apastron (dotted lines) are calculated for initial 
Lorentz factor T = 10 3 , 10 4 , 10 5 , 10 6 and 10 7 (indicated at the 
curves). The luminosity of the companion star was assumed 
^star = 2.2 ■ 10 38 erg/s. The experimental points correspond 
to the average energy spectrum measured by HESS within 
several weeks around the periastron, note that the last mea- 
surement was multiplied by a factor of 10. The upper limits 
in energy r ange 30MeV-3GeV correspond to the EGRET mea- 
surements (Ta vani et al,|[T996l) . The differential flux sensitivity 
of GLAST for a point source is also shown. The corresponding 
curve represents sensitivity for one-year all-sky survey taken from 
http: / /www-glast. slac.stanford.edu/software/IS/glast_lat_performance. htm 
However, since gamma-ray fluxes can be effectively observed only 
at the epochs not far from the periastron, the typical available 
observation time by GLAST would be limited by At ^ 3 weeks. 
This implies that the minimum detectable gamma-ray flux of 
GLAST shown in the figure should be increased by a factor 
between (lyear/At) 1 / 2 __ 4 and 16 depending whether the 
sensitivity is determined by the background or by the gamma-ray 
photon statistics. 



